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We present a measurement of the absolute branching fraction for D^ -^ 

1^ • K~TT^ using the reconstruction of the decay chain B -^ D*^Xi~r', D*^ — > 

Cn ' ZJ^TT"*" where only the lepton and the low-momentum pion from the D*^ are 

detected. With data collected by the CLEO II detector at the Cornell Electron 

CN ! Storage Ring, we have determined Br{D^ -^ K^-k^) = [3.81 ± 0.15{stat.) ± 

0.16{syst.)]%. 
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As most of the published branching fractions of D^, D^ and D^ mesons are normahzed 
to the D*^ -^ K^Ti^ |jl|] decay mode, then the value of Br{D^ — > K~tt^) directly affects 
many topics in heavy flavor physics. Some examples include charm counting in B meson 
decays where about 90% of the total charm yield is calibrated by Br{D^ -^ K~7i~^) 0, the 
determination of Br{Z^ -^ cc), and the investigation of any exclusive decay mode of the B 
meson which contains D^, D^ or D^ in the final state. 

In order to measure the absolute branching fraction for D^ — * K^tt^ decay, one needs to 
find the number of D°'s without reconstructing a particular D^ decay mode. In this Letter 
we present a measurement of the absolute D^ — > K~tt~^ branching fraction, developing 
the method first used by the ARGUS Collaboration ||^. The inclusive number of Z^^'s is 
determined by partial reconstruction of the decay chain B —>■ D*^i~iJ, D*^ -^ D^n^, 
where only the lepton and the slow pion from the D*~^, hereafter denoted as tCs, are detected. 
The systematic errors involved are largely different from those of other recent measurements 
10-^, where slow pions within jets were used to tag the decay D*^ -^ D^ti^ . 

We have used 3.1 /6~^ of data collected on the T(4S') resonance by the CLEO II detector 
0. The data set corresponds to 3.3 x 10^ BB events. In order to suppress non-BB (contin- 
uum) background we required the ratio of the Fox- Wolfram moments H2/ Hq [Q to be less 
than 0.4. The remaining contribution from continuum events was estimated using 1.6 fh~^ 
of data collected just below the BB threshold. In the following this continuum subtraction 
is implicit. 

We required lepton candidates to have a momentum between 1.4 GeV/c and 2.5 GeV/c 
and to be in the barrel region of the detector. Muon candidates were required to penetrate 
an iron absorber to a depth of at least 5 nuclear interaction lengths. Electrons were identified 
through a comparison of the energy deposited in the electromagnetic calorimeter with the 
momentum measured in the drift chambers and by specific ionization energy loss (dE/dx) 
measurements. We required that the vr^ candidate have a momentum lower than 190 MeV/c, 
which is slightly below the upper kinematic limit for pions from D*^ m B ^ D*'^i~u decays. 

The partial reconstruction of the decay B -^ D*~^i~iJ exploits the extremely low energy 
release in the decay D*~^ — >■ D^tt^ . The pion is almost at rest in the D*^ frame, and its 
velocity vector in the lab frame is approximately equal to that of the D*^ . Our main signal 
mode is B -^ D*^i^jj, for which the missing mass squared is calculated as 

MM^ = {Eb -Ee- En^+f - \Pb -Pi- Pd*+\^- (1) 

The energy of the B meson is precisely the beam energy. We do not know the direction of 
motion of the B, but the B momentum is sufficiently small ( ~ 300 MeV/c) compared to 
the typical values of \Pi\ and |-Pd*+| that we can set Pb = 0. We approximated the direction 
of motion of the D*^ by the direction of motion of the vr^. If the tt^ were exactly at rest in 
the D*^ frame, the D*^ energy would be given by E''£i+ = {E^^^ / E'j;^-) ■ Mj:,*+. In order to 
correct for the non-zero momentum of the vr^ in the D*^ frame, we used a parameterization 
obtained from Monte Carlo to estimate Ej:,*+ as a function the vr^ momentum |^. 

The resulting MM"^ distribution is shown in Figure |l|(a). The events with the lepton 
and slow pion coming from B -^ D*~^i~i', D*^ -^ D'^'^t produce a prominent peak at 
MM"^ ^ 0. However, the decays B -^ D*^X£~u, D*^ -^ -D°7r+ also contribute to this peak. 
We have considered these decay modes to be signal because they produce true D*~^ — > D^nf. 



More specifically, we allowed the D*^ to come from B -^ D*^n7T£~u decays, where D*^n7i 
may or may not form a resonance. We also allowed the lepton to come from r in the decays 
B -^ D*^T~i' or from D in the decays B -^ D*^DX, where D represents D , D~ or D~ . 
Our analysis is therefore not dependent on the branching fractions assumed in the Monte 
Carlo for the poorly measured B -^ D*'^n7ri~i' and B — > D*^DX decays, because these 
decays were considered to be signal. 

A Monte-Carlo simulation of the BB events was used to determine the background 
shape. We normalized the background shape to the data distribution in the sideband region 
[MM"^ < —5 GeV^/c^). After the background subtraction, the number of events in the signal 
region (defined as MM^ > -2 GeVVc^ ) was found to be A^*"'^' = 44, 504 ± 360 (stat.). In 
this way we have extracted the number of i? ^ D*^Xi~i' events in which D*^ -^ -D°7r+. 

We have thus obtained a sample of D*^ -^ D^n^ decays without reconstructing a partic- 
ular D^ decay mode. Next we need to determine how many Z^^'s from these D*~^ — *■ D^7r~^ 
events decay to K~7t~^. For every i~7r^ pair for which the value of MM^ was within the sig- 
nal region {MM"^ > —2 GeV^/c^) we searched for a K^n^ pair, assigning the kaon mass to 
the track of the opposite charge with respect to tt^, and requiring \M{K~'k'^) —M{D^)\ < 35 
MeV/c^, which corresponds to a 3.5 cr cut. The K^tt^ pair was combined with the n^ and 
the mass difference AM = M^K^ti^tt^) — M{K~'tt^) was formed. The resulting AM dis- 
tribution is shown in Figure §. The prominent peak at AM = M{D*^) — M{D^) ^ 145.4 
MeV/c^ is produced by D*~^ — > D^n^, D^ -^ K^tt^ decays. We normalized the background 
shape obtained from the Monte-Carlo simulation to the data distribution in the sideband 
region (155 MeV/c^ < AM <180 MeV/c^). True D*+ -^ D^tt+, D° -^ K-n+ decays where 
the D*~^ does not come from a signal decay chain were considered to be background. After 
the background subtraction we counted the number of events in the signal region, defined 
as 141.50 MeV/c2 < AM <149.75 MeV/c^. The number of decays D*+ -> D^tt+ with 
D^ -^ A'^7r+, denoted as N^^^\ was found to be 1165 ± 45 (stat.). 

To extract Br{D^ —>■ K^-k^) we need to correct the ratio jsf'^^'^^ j fsfi'^'^^ for the track 
reconstruction and acceptance efficiencies : 

Tvrexcl -| 

BriD"" -^ K-n+) = —— ■ -. (2) 

We obtained e using a GEANT-based Monte-Carlo simulation 0] of the CLEO II detec- 
tor. To a good approximation the lepton and slow pion reconstruction efficiencies cancel 
in the ratio when we calculate e. Therefore e mainly includes reconstruction and selection 
efficiencies for K~ and tt"*" tracks and acceptance efficiencies for the M{Kt[) and AM signal 
regions. However, the cancellation of the lepton and slow pion reconstruction efficiencies is 
not exact because the average charged track multiplicity for D° decays is higher than that 
for D^ —>■ K^TT^ mode and it is more difficult to reconstruct a track in a higher multiplic- 
ity environment. We found that this effect changes e by 3.7% of itself. In order to take 
this into account, we calculated e by selecting signal events from the Monte-Carlo simula- 
tion of BB events, and comparing the value of N'^^ / N^^f, to the branching ratio that was 
used in the Monte Carlo. Note that in this procedure A"|,j^ corresponds to the number of 
B -^ D*~^X£^u, D*^ -^ D^TT^ events where Z)°'s were allowed to decay generically, not 
forced to decay into K^-n^ . We obtained e = [68.6 ± 2.1{syst.)]%, and using this value of e 
together with Eqn. ^ we found 



Br{D° -^ K'Ti^) = [3.81 ± 0.15{stat.) ± 0.16{syst. 



0. 



The total systematic error was obtained by summing in quadrature the errors given in 
Table I. We will now discuss the systematic uncertainties dividing the possible sources into 
three categories: (i) determination of A^*"'^' using the MM"^ distribution, (ii) determination 
of N'^^'^^ using the AM distribution, (iii) efficiency extraction from Monte Carlo. 

(i) First, to see how well the Monte Carlo can simulate the background shape for the 
MM"^ distribution, we looked at the MM"^ distribution for the wrong-sign (i.e. same sign) 
ills pairs (Figure |l](b) ). We normalized the Monte-Carlo shape to data distribution in the 
sideband region {MM"^ < — 5 GeV^/c^), as we did for the right-sign iiTg pairs, and compared 
the Monte-Carlo prediction with data in the signal region {MM"^ > —2 GeV^/c'^). We found 
excellent agreement within the statistical precision of 0.8% of the signal region population. 
We include this 0.8% as a part of the systematic error. This result is encouraging, but differ- 
ent physics can contribute to the distributions for wrong-sign and right-sign background £ vr^ 
pairs. Using Monte Carlo, we performed a thorough study comparing the MM"^ distributions 
for the various physical processes producing the wrong-sign or the right-sign background ^vr^ 
pairs. 

We have found that the most dangerous source of background which peaks in the signal 
region of MM"^ distribution is the decay chain B — > DXi~iJ, D -^ (something heavy) -|- vr"*", 
where the Tr"*" is moving slowly in the D rest frame and mimics the pion from D*^ -^ D^T^t 
decay. These decays do not contribute to the AM peak and thus can reduce the measured 
D^ — ^ K~Ti^ branching fraction. To estimate the systematic error due to this background 
we identified three such low Q-value decay modes in our Monte Carlo: D^ — ^ K (892)u;7r+, 
D^ — > K p'^TT'^, and D^ — > K p°7r+. Monte Carlo predicts that the events with the pion 
coming from one of these modes account for 0.7% of the events under the MM'^ peak with 
respect to the number of events in the signal peak. We have exploited the difference in 
the MM"^ distribution shapes for this background and the signal and fit the whole MM^ 
data distribution with three histograms obtained from Monte Carlo: signal, the contribution 
from the decay chain B —>■ D^X£~u where D~^ —>■ K [u: or p)7T~^ , and the rest of background. 
The fit showed that the contribution from these modes is consistent with the Monte-Carlo 
prediction. However we should keep in mind that the decay modes we are considering here are 
poorly measured and that there could be other similar low Q-value decays that have not yet 
been observed. In order to be conservative, we varied the contribution from B — > D~^Xi~i?, 
D^ —>■ K {uj or p)n~^ in the Monte-Carlo background shape by the fit error and obtained a 
2.3% variation in final result, which we took as the systematic error due to this background. 
This is the largest single source of systematic uncertainty in the analysis. 

Another source of background which peaks in the signal region of the MM'^ distribution 
results when the slow pion from a signal decay chain decays in fiight to a muon, and we 
identify this muon as the slow pion. Monte Carlo predicts the magnitude of background 
from this source in the MM"^ peak region to be 2.5% of the signal. Even though this is the 
largest source of background which peaks in the signal region it does not significantly bias 
the Br{D^ —* K~7c~^) measurement because this background produces smeared peaks in the 
signal regions of both the MM"^ and the AM distributions. We varied the Monte-Carlo 
prediction for this background by 30% of itself and obtained 0.3% variation in final result, 
which we took as the systematic error. 
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Another background which peaks in the MM^ signal region results when we identify 
as a 7r+ a positron from 7r° -^ 7e+e^ or 7 conversion in the decay chain B — > D*Xi~u, 
D* — ^ Dtt^, D'-f. Monte Carlo predicts the magnitude of background from this source in the 
MM"^ peak region to be 0.7% of the signal. We varied the Monte Carlo prediction for this 
background by 30% of itself and obtained 0.4% variation in final result, which we took as 
the systematic error. 

Combining the errors described above in (i) we estimated the systematic error due to 
background subtraction in the the MM'^ distribution to be 2.5%. We have also studied 
the possible systematic errors due to the cut on slow pion momentum, fitting and yield 
determination in MM"^ distribution, and fake leptons. The results of these studies are given 
in Table I. 

(ii) We have studied the systematic error due to the background subtraction in the AM 
distribution. We included true D*^ — ^ D^^^t ■• -^° ~^ K~t[^ decays where the D*^ does not 
come from a signal decay chain in the definition of background. The main source of this 
background is D*'^i~ pairs for which the D*^ comes from one B , and the lepton is the 
primary lepton from another B . This background is suppressed because it occurs only due 
to B^ — B mixing. A less significant source is D*^i^ pairs for which the D*^ comes from one 
B or B^ and the lepton is a secondary lepton from the D from the other B^ or B~^. This 
background is suppressed by the lepton momentum requirement which predominantly selects 
primary leptons from B decays. Neither of these background components contribute to the 
peak at MM'^ ^ because the lepton and slow pion come from different S's. We varied the 
Monte Carlo prediction for these backgrounds by 20% (based on the conservative estimate of 
the uncertainties in the inclusive D*^ and lepton yields, the B^ — B mixing parameter, and 
the dependence of MM"^ distribution shape on the D*^ momentum spectrum), and obtained 
0.6% variation in final result, which we took as the systematic error. 

The rest of the background in the AM distribution is combinatoric. To estimate the 
systematic error due to the Monte Carlo simulation of this background we substituted the 
combinatoric part of the Monte Carlo background shape by an analytic threshold function 
and obtained the 0.9% shift in the final result, which we took as the systematic error. 

Combining the errors described above in (ii) we estimated the systematic error due to 
background subtraction in the the AM distribution to be 1.1%. We have also studied the 
possible systematic error due to the fitting and yield determination in the AM distribution, 
and the result of this study is given in Table I. 

(iii) A study has been performed to estimate the systematic error due to the extraction 
of the reconstruction efficiency for K~ and 7r+ tracks from Monte Carlo. We assigned a 
2% error to the final result (1% per track). As was mentioned earlier, the lepton and slow 
pion reconstruction efficiencies do not cancel out exactly due to the difference in charged 
multiplicity between the cases D^ — * K^-k^ and D^ -^ all. To estimate the systematic error 
due to this effect we extracted the efficiency from Monte Carlo forcing D^ -^ K'-k^ when 
we determine N'^mc- -^^ ^ systematic error we took 30% of the shift in the efficiency obtained 
using this method and the method actually employed in the analysis. We have also studied 
the possible systematic error due to the choice of the signal region in the AM distribution, 
and the result of this study is given in Table I. 

The systematic errors due to the limited Monte Carlo statistics and the continuum sub- 



traction are also given in Table I. 



Quantity Possible source of systematic error 


Estimate of Error 
(% of final result) 


N'^"-'^^ Background subtraction in MM^ distribution 

Slow pion momentum cut (affects MM^ background shape) 
Fitting and yield determination 
Fake leptons 


2.5% 
1.0% 
0.6% 
0.2% 


jyexd Background subtraction in AM distribution 
Fitting and yield determination 


1.1% 
0.3% 


e K^-K^ reconstruction efficiency 

Choice of signal region in AM distribution 

Non-exact cancellation of i and tTs reconstruction efficiencies 


2.0% 
1.6% 
1.1% 


Monte Carlo statistics 


1.4% 


Continuum subtraction 


0.1% 


Total 


4.3% 



TABLE I. Systematic error summary table. 



In conclusion, we have measured the absolute branching fraction for D^ ^ K n^ decay 
using a 5 -^ D*+Xtp tag. We have found Br{D^ -^ K'tt+) = [3.81 ± 0.15{stat.) ± 
0.16{syst.)]% |TT]]. Our result is consistent with a recent measurement by ALEPH of (3.82 ± 
0.09 ± 0.11)% i, g two measurements by ARGUS of (3.41 ± 0.12 ± 0.28)% g and of 
(4.5 ± 0.6 ± 0.4)% §, and two measurements by CLEO of (3.91 ± 0.08 ± 0.17)% § and of 
(3.69±0.11±0.16)% ||12|. Taking into account correlations, we combined our result with the 
other two CLEO measurements and found a new CLEO average value for Br{D^ -^ K^tt^) 
to be [3.82 ± 0.07{stat.) ± 0. 12 {syst.)]%. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, the A. P. Sloan Foundation, and the Swiss National Science Foundation. 



^We took the value before correction for the final state radiation from the K and n daughters in 
the D^ decay. 
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FIG. 1. The missing mass squared (MM'^) distribution for the right-sign (a) and wrong-sign 
(b) iiTs pairs. The Monte Carfo background shape has been normahzed to the data distribution in 
the sideband region indicated by the dashed hne (MM'^ < —5 GeV^/c^). 
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